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Abstract 

Theory of the optical parametric amplification at high-frequency pumping in crystals 
with a regular space modulation of the sign of nonlinear coupling coefficient of interacting 
waves is developed. By applying the matrix method, the theory is based on a step-by-step 
approach. It is shown that, in the case where the pumping intensity is less than some 
critical value, the spatial dynamics of the signal intensity inside a separate layer with the 
constant nonlinear coefficient has an oscillatory behavior and the change of the signal 
intensity from layer to layer is defined, in general, by the power function. The same law 
is valid for the change of variance of signal's quadrature components. At large number 
of layers, these dependences can be reduced to the well-known ones for homogeneous 
nonlinear optical crystals. 

KEY WORDS: parametric interaction, periodically poled nonlinear crystals, quasi- 
phase matching, step-by-step approach, quadrature component 

1 Introduction 

Quasi-phase-matched interactions of light waves as well as phase-matched interactions 
provide the realization of an effective energy exchange between interacting waves (see, 
for example, the reviews [H and the possibility to obtain nonclassical light E]. 



Quasi-phase-matched interactions are implemented in crystals with the periodic modu- 
lation of nonlinear susceptibility. These are crystals with the regular domain structure 
or periodically-poled nonlinear crystals (PPNC). In such crystals, the mismatch of wave 
vectors Ak of interacting waves is compensated by an inverse-lattice vector k, of nonlinear 
susceptibility's modulation 

K = 271171/ A, 

where A is the modulation period and m is the quasi-phase-matching order. 

The conventional approach for implementing parametric interactions in PPNC consists 
in a choice of the domain-structure period A, in view of the expression 

A = 2n/\Ak\, 

where Ak is the phase mismatch between interacting waves. The period A is usually 
determined from the analysis of the second-harmonic generation |2] . At the coherence 
length /c = A/2, the second-harmonic intensity reaches the maximum value and the phase 
(related to the phase mismatch and equal to vr) compensates the change of the nonlinearity 
sign in the neighboring layer (domain) due to the optical-axis inversion. 

While nonlinear optical interactions are calculated analytically, PPNC is replaced by a 
homogeneous crystal with an effective nonlinear coefficient jSes distinguished from a value 
(3 for the homogeneous crystal by the factor 2/7rm, so that 

f3es = —f3. 

However, such approach to the optical-parametric-amplification process is not always 
valued. As it is known, in a parametric process the transition to a steady state of the 
optimum phase relation essentially depends on the pumping intensity |S1IS1- Recently, the 
spatial dynamics of intensity and phase of the parametrically amplified wave in PPNC 
was investigated in detail |7]. In particular, it was shown that at the transition from 
layer to layer the sign of derivative in the phase relation is changed. The change of the 
signal wave intensity in PPNC also differs from the homogeneous crystal case. In the 
present paper, the simple "step-by-step" approach to the analysis of quasi-phase-matched 
optical parametric amplification is presented. The knowledge of transmission factors of the 
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layers is the basis of this approach. For the definition of transmission factors, the matrix 
approach elaborated in |8,; is used. For an arbitrary number of layers with thickness equal 
to the coherence length, expressions for the signal intensity and variance of quadrature 
components for an initial random signal phase are derived. 
The paper is organized as follows. 

In Sec. 2 the equations for degenerate optical parametric amplification in PPNC are 
presented in generic form using undepleted pump approximation. Section 3 contains the 
solution of differential equation for the optical parametric amplification in PPNC in ma- 
trix form. In Sec. 4 the solution obtained is analyzed for the case satisfying the condition 
of quasi-phase-matched interactions. In Sec. 5 the possibility to generate the quadrature- 
squeezed light is considered when classical fiuctuations are suppressed. The results ob- 
tained are summarized in Sec. 6. 

2 Basic Equations 

The optical parametric amplification process in PPNC of a wave with frequency uo in the 
pumping-wave field with frequency 2uj is determined by the following equations: 



where A and Ap are the complex amplitudes of the signal and pumping waves, respec- 
tively, (3 is the modulus of nonlinear wave coupling coefficient, Ak = k — 2kp is the phase 
mismatch, and g{z) is the periodic function equal either to -|-1 or to —1 at the thickness 
of a separate layer, being dependent on the nonlinear-coefficient sign. In the general case, 
it is impossible to solve the set of equations ((H) analytically. 

Within the framework of undepleted pump approximation, the process under consid- 
eration is described by the equation 



where Ap is the constant value. Equation Q makes it possible to get an analytical solution 
on each layer for given g{z). However, the solution for arbitrary number of layers, cannot 




■if3g{z)ApA* exp{—iAkz), 
-iPg{z)A'^ exi){iAkz), 





(2) 
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be cast into a suitable for the application form. In this connection, to solve Eq. Q, we use 
the matrix method, which allows one to obtain the solution in a convenient form. Firstly 
we rewrite Eq. Q introducing a reduced length ( = z/Lj^i. 

dA 

— = -z(7(C)A*exp(-z5C), 

(3) 

A{C = 0) = Ao, 

where Lni = l/Z^l^pl is the so-called nonlinear length, 5 = AkL^\ is the normalized phase 
mismatch, and Aq is the signal amplitude at the PPNC input. 

3 Matrix Differential Equation and Its Discrete So- 
lution 

Equations (jH} by means of the substitution 

A = Bexp{-i6C/2) (4) 

are reduced to the following system of equations: 

dB i5B . 



^ = -— + W(C)B. 
which can be presented in matrix form 

with the initial values C = Cq. In Eq. (6) Cq is the value of the matrix C at C = 
and D is the matrix determined by the expression 

^(0 = 



t6/2 -tgiC) 
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The function g{() takes the value either +1 or —1. One introduces matrices D± associated 
with these values 
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D± = 7^\ \ (9) 

^ » ±2 -5 

Let PPNC have layers with the length A/2. It is just the case used for repolarization 
of a nonlinear crystal. Let g{() = 1 in the first layer. Then the solution of Eq. (jH)) reads 

C(C) = exp{D^OCo, (10) 

where the exponential function of matrix -D+C is nothing else as the transmission factor 
for the layer and it is introduced as follows: 

-p(«+o = E^- (") 

k=Q 

In accordance with (fTUI) . at the output of the first layer, we have 

C(A/2) = exp(D+A/2)Co. (12) 

In order to obtain the signal at the output of the second layer, we need to suppose D = 
and take into account the fact that the output of the first layer is simultaneously the input 
of the second layer. As a result, the signal at the output of the second layer is determined 
by the expression 

C(A) = exp(D_A/2)C(A/2) = exp(D_A/2) exp(D+A/2)Co. (13) 
For layer pairs, Eq. (13) reads as follows: 

r 1^ 

CiNA)= exp(D_A/2) exp(D+A/2) Co = Q^Co, (14) 
where the matrix Q is 

Q = exp(D_A/2) exp(L'+A/2). 

Thus, the matrix defines the signal transformation from the PPNC input to its output. 
The solution constructed by such a way is based on the assumption that PPNC has a 
discrete structure. Due to this, we call such approach the step-by-step approach. 
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4 Quasi-Phase-Matched Interaction 

Now we calculate the matrix Q in the case where the quasi-phase- matching condition is 
realized and the modulation period A = 2/c- The exponential matrix can be presented in 
a compact form 

(15) 

±2is c — ids I 

where functions c(A) and s(A) read 

c(A) = cos(/iA/4), s(A) = /i"^ sin(/iA/4), (16) 

with 



In these relations, Ip = |ApP is the pumping intensity, = (AA;/2/3)^ is the so-called 
critical intensity of the pumping wave or the parametric-trapping intensity. To use PPNC 
for the optical parametric amplification makes sense at Ip < lev In this case, the signal 
wave inside a separate layer oscillates. By multiplying the exponential matrices D±, the 
matrix Q can obtained in the form 

I -As^ + ic + is5f As^5 \ 

Q=\ • (17) 

y As^5 -As^ + (c - is5f j 

Since A = 21^ = 2n/\6\, the functions (fTBj) are 

c(g=cos(|0, .(y =/.-isin(^0. (18) 

At \6\ ^ 2 (Jcr ^ Ip), fJ'/\S\ ~ 1, therefore, we arrive at 

1 / -4-52 4151 

Q^--2\ \ (19) 

\^ A\S\ -4-52 

Raising the matrix (fT^ to A^th power, we obtain 

1. «^(^|. (20) 



Taking into account expressions (jZj) and p4j] . we arrive at the following result for B: 

(a-^ + a^)5o + {a-'' - . (21) 
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Now in view of the relations 

Aq = aoe'"^" and Iq = |Ao|^ 



we obtain the following expression for the signal wave intensity 1^ = \B\'^ = \A\'^ at the 
PPNC output: 



IN 



-2N 2 , 2N ■ 2 

a cos (fo + a sm ipo 



Jo. (22) 

One can see that it is just the well-known property of phase sensitivity of degenerate 
parametric amplification. From Eq. ()22|) follows that the largest amplification takes place 
at the signal phase ipo = ±(7r/2). The largest signal suppression is achieved at the signal 
phase ipo = ivr. For number of layer pairs ^ 1, at ^ 2 expression (j^^ can be 
transformed into the following one: 



cosh(rL) - cos(2(/?o) sinh(rL) Iq, (23) 



where F = (4/7r)/5|74p| is the parametric amplification increment and L = 2Nlc is the total 
length of PPNC. It is worthy to compare expression (f^!?|) with the analogous one for the 
case of a homogeneous crystal under the phase-matching condition. In the latter case, 
the parametric amplification increment is equal to Fq = 2/?o|^p| (see, for example, [H]), 
where (3o is the nonlinear wave coupling coefficient for the homogeneous crystal. The ratio 
F/Fq = {2 / it) f3 / (3q shows the difference between the optical parametric amplification in 
PPNC and in the homogeneous crystal. In the case f3 = (3q, we obtain the well-known 
result F/Fo = 2/7r. 

Thus, at the quasi-phase-matched-wave interaction, the intensity of parametrically 
amplified signal wave under the conditions mentioned above varies in the same way as in 
the homogeneous crystal with the effective nonlinear coefficient f3es = (2/7r)/?. It should 
be noted that the same value (3es for PPNC can be obtained, if truncated equation is 
averaged over the modulation period of nonlinear- wave-coupling coefficient. However, in 
the latter case the conditions of applicability of the equations obtained are not evident. 



Using formula (j221) we have calculated curves presented in Figs. 1 and 2 for ip^ = 
tt/2. Figure 1 shows the dependence of the signal intensity on the interaction length at 
the quasi-phase-matched parametric amplification in PPNC Figure 2 demonstrates the 
dependence of the parametric amplification increment on the interaction length at different 
relations between the pumping intensity and its critical value. From these figures, one 
can see that, at the same pumping intensity and given interaction length, the signal wave 
intensity increases when the pumping intensity approaches its critical value. Therefore, the 
replacement of PPNC by the homogeneous crystal with the effective nonlinear coefficient 
provides less value for the signal intensity. This difference grows up when the pumping 
intensity Ip approaches to /„. Comparison of the conversion efficiency for the parametric 
amplification in PPNC and in the homogeneous crystal with nonlinear coefficient Pes is 
done in Fig. 3. Note that at Ip ^ /„ we have used general formula (fT7|). 

5 Quadrature Components 

We turn now to the analysis of quadrature components X, Y of the signal wave in PPNC: 
X{z) = \ [A{z) + A\z)] , Y{z) = 1 [A{z) - A*{z)] . (24) 
According to Eq. (PT|) . one has for X-quadrature 

X{z) = ^ [{a-^ + a^) cos 0_ + (a^^ - a^) cos 0+] , (25) 

where the phase (l)± = (po ± 6^/2. 

To obtain a signal with suppressed fluctuations of one of the quadrature components, 
the random signal with a uniform distribution function for the initial phase (fo is worthy 
to consider 

= 77-' -VT < v^o < TT. (26) 

ZTT 

Then the mean value {X{z)) = and the variance of quadrature component 
{X'^{z)) = X'^{z) — {X{z))'^ is given by the expression 

{X'iz)) = ^-al [a^^ + (a-^^ - e^^) cos^(5C/2)] . (27) 
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At S( = AkL = 2Tin (n = ±1, ±2, . . .), X-quadrature fluctuations are reduced, namely, 

{X\z))^,, = (28) 

In the case 5C, = 'n:{l + 2n), the variance reaches the maximum value: 

(X2(^))max = la'y, (29) 

i.e., fluctuations of quadrature components are amplifled. 

Thus, in PPNC the spatial changing of the variance of signal's quadrature components 
is determined by the power functions as well as the signal intensity. However, if number 
of layer pairs N ^ 1, expressions ()27 |) -()29 |) can be transformed into the forms similar to 
the ones for the homogeneous crystal: 

(X2(^))min/max = {X\0))e^^'^. (30) 

It is not difficult to draw a conclusion, in view of formulas p8|l and (j30|l . that suppression 
of fluctuations in PPNC occurs more effectively than in the homogeneous optical crystal 
with the effective nonlinear coefficient. 

6 Conclusions 

The application of step-by-step matrix approach allows one to calculate the optical para- 
metric amplification in periodically-poled nonlinear crystals without analysis of the spatial 
dynamics inside a separate layer. The method can be applied to analysis of the processes 
in PPNC at an arbitrary period of the crystal repolarization and arbitrary relation of the 
critical intensity to the pumping intensity. In other words, the method is not restricted 
to the conditions Ip <^ and N ^ 1 which are necessary for the correct analysis in 
the traditional-approach case where PPNC is replaced by an homogeneous crystal with 
effective nonlinear coefficient. This circumstance allows one to use this method for both 
analytical and numerical [Hj calculations. 

The spatial dynamics of the signal intensity and quadrature components was analyzed 
versus the ratio Ip/Ia and the interaction length. It was shown that at the pumping 
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intensity less than the critical value, the spatial distribution of the signal intensity inside 
a separate layer exhibits an oscillatory character and changing of the signal intensity 
and variance of quadrature components from layer to layer follow the power function 
dependences. In the case of larger number of layers and Ip <^ /cr, these dependences can 
be reduced to well-known formulas for homogeneous nonlinear crystals. Notice that the 
requirement Ip > Icr corresponds to the parametric trapping in the homogeneous crystal. 

It should be emphasized that our analyze has shown that the quasi-phase-matched 
optical parametric process in PPNC occurs, generally speaking, more effectively than the 
parametric process in the homogeneous crystal with the nonlinear coupling coefficient 
equal to the effective nonlinear coefficient for PPNC. The generalization of the approach 
presented here onto the case of quasi-phase-matched consecutive wave interactions and 
nonstationary parametric interactions is in progress. 
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Figure 1: Coefficient of the parametric signal amplification in PPNC as a function of the in- 
teraction length for the phase ipo = 7r/2 and different ratios of the pumping intensity Ip to the 
critical intensity /„: Ip = Icr (1), Ip < Icr (2), and Ip <C Icr (3). Curve 2 corresponds to formula 
H22|) and curve 3 corresponds to formula ()23() . 
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Figure 2: Relative increment of the parametric amplification in PPNC as a function of the 
interaction length and ratio I^r/Ip- I and /he are the signal intensities for the case of PPNC and 
homogeneous nonlinear crystal, respectively. 
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Figure 3: Ratio of the parametrically-amplified signal intensity / in a LiNbOs PPNC to the 
intensity Ihc in a LiNbOa homogeneous crystal with corresponding effective nonlinear coefficient 
/3efr versus 7 = fiyJTp for the eee-interaction type and A = 0.5 ;um, A = 20 fim, N = 3, and 
7cr = /?\/^- 
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